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We report the detection of a substellar companion orbiting the intermediate- 
mass giant star 11 Com (G8 III). Precise Doppler measurements of the star from 
Xinglong station and Okayama Astrophysical Observatory (0A0) revealed Kep- 
lerian velocity variations with an orbital period of 326.03±0.32 days, a semiampli- 
tude of 302.8±2.6 m s _1 , and an eccentricity of 0.231±0.005. Adopting a stellar 
mass of 2.7±0.3 M Q , the minimum mass of the companion is 19.4±1.5 Mj, well 
above the deuterium burning limit, and the semimajor axis is 1.29±0.05 AU. This 
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is the first result from the joint planet search program between China and Japan 
aiming at revealing statistics of substellar companions around intermediate-mass 
giants. 11 Com b emerged from 300 targets of the planet search program at OAO. 
The current detection rate of a brown dwarf candidate seems to be comparable 
to that around solar-type stars within orbital separations of ~3 AU. 

Subject headings: stars: individual: 11 Com (HD 107383) — stars: low- mass, 
brown dwarfs — planetary systems — techniques: radial velocities 

1. Introduction 

A "brown dwarf desert" is widely known as a deficit in the frequency of brown dwarf 
companions (13 — 80Mj) to solar-type stars (e.g. Marcy & Butler 2000). The mass distri- 
bution of companions rises steeply toward planets at one extreme and stellar masses at the 
other, with a deficit of brown dwarfs in between (e.g. Grether & Lineweaver 2006). The 
presence of such a desert within orbital separation of ~3 AU was revealed and confirmed by 
the precise radial velocity surveys during the past decade (Butler et al. 2006 and references 
therein), and recent coronagraphic surveys showed that the desert extends to ~1000 AU 
(McCarthy & Zuckerman 2004). This bimodal mass distribution of the low-mass compan- 
ions is considered to be the outcome of distinct formation mechanisms for planets and stellar 
companions to solar-type stars. 

On the other hand, it has remained to be cleared whether the brown dwarf desert also 
exists around massive stars. Few absorption lines in the spectra of massive early-type dwarfs, 
which are often rotationally broadened, make it difficult to detect low-mass companions by 
precise radial velocity measurements, and their high luminosity prevents us from directly 
imaging close and faint companions. Therefore, massive stars had not been major targets of 
planet/brown dwarf searches compared with lower-mass stars. 

However, the existence of substellar companions to intermediate- mass stars (> 1.6M Q ) 
has been gradually uncovered mainly through the precise Doppler surveys of G, K giants 
(e.g. Setiawan et al. 2005; Sato et al. 2007; Hatzes et al. 2005; Johnson et al. 2007; 
Lovis & Mayor 2007; Frink et al. 2002; Niedzielski et al. 2007). These are massive stars in 
evolved stages and have many sharp absorption lines in their spectra appropriate for radial 
velocity measurements. The minimum masses of the discovered companions range from 0.6 
to 19.8 Mj, and the semimajor axes are from 0.8 to 2.4 AU. The two of the most massive 
ones, NGC4349 Nol27 b (19.8 My Lovis & Mayor 2007) and HD 13189 b (14 My Hatzes 



- 3- 



et al. 2005) ft are falling in the brown dwarf regime. Galland et al. (2005) have tried a 
precise Doppler survey of low-mass companions to A-F dwarfs and found one brown dwarf 
candidate with the minimum mass of 25Mj in a 25-day orbit around the A9 V star HD 
180777. Kouwenhoven et al. (2007) found two distant (520 and 1500 AU) brown dwarf 
candidates to a B star in the Sco OB2 association by direct imaging with the near-infrared 
camera and the adaptive optics system. They derived a frequency of brown dwarf companions 
(> 30Mj) to 0.5±0.5% in the separation range 130-520 AU for late-B and A type stars in 
the association, which is comparable to the frequency of 0.7±0.7% derived by McCarthy & 
Zuckerman (2004) for F-M stars in the separation range 120-1200 AU. 

We here report the detection of a substellar companion to the intermediate-mass G-type 
giant 11 Com (HD 107383) from Doppler planet search programs at the National Astronomi- 
cal Observatories (Xinglong, China) and Okayama Astrophysical Observatory (OAO, Japan). 
This is the third one with the minimum mass falling in the brown dwarf regime discovered 
around intermediate-mass giants. 

2. Observations 

2.1. OAO Observations 

The Okayama Planet Search Program started in 2001 using a 1.88 m telescope and the 
High Dispersion Echelle Spectrograph (HIDES; Izumiura 1999) at OAO. It aims to detect 
planets around intermediate-mass G-type (and early K-type) giants (Sato et al. 2005) and 
about 300 targets are now under survey. For precise radial velocity measurements, we set 
a wavelength range to 5000-6100A and a slit width to 200 /im (0.76") giving a spectral 
resolution (A/AA) of 67000, and use an iodine absorption cell (I2 cell; Kambe et al. 2002) for 
precise wavelength calibration. Our modeling technique of an I 2 -superposed stellar spectrum 
(star+I 2 ) is detailed in Sato et al. (2002), which is based on the method by Butler et al. 
(1996), giving a Doppler precision of about 6 m s _1 over a time span of 6 years. The stellar 
template used for radial velocity analysis is extracted from several star+12 spectra with use 
of the method described in Sato et al. (2002). The reduction of echelle data is performed 
using the IRAF I software package in the standard manner. For abundance analysis of 11 



1 HD 13189 b has uncertainty between 8 to 20 Mj depending on the uncertainty of the host star's mass 
(2-6 M Q ) (Hatzes et al. 2005) 

2 IRAF is distributed by the National Optical Astronomy Observatories, which is operated by the Associ- 
ation of Universities for Research in Astronomy, Inc. under cooperative agreement with the National Science 
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Com, we took pure (I 2 -free) stellar spectra of wavebands of 5000-6100 and 6000-7100 A with 
the same wavelength resolution as that for radial velocity measurements. 

2.2. Xinglong Observations 

The Xinglong Planet Search Program started in 2005 within a framework of international 
collaboration between China and Japan aiming to extend the planet search program ongoing 
at OAO. At the National Astronomical Observatories (Xinglong, China), we have mainly 
monitored about 100 G-type giants with a visual magnitude of V ~ 6, which have not been 
observed at OAO. 

Before starting the planet search program, we installed an I 2 cell to the Coude Echelle 
Spectrograph (CES; Zhao & Li 2001), which attached to 2.16 m telescope at Xinglong 
Station in 2004 August. The cell is a copy of those for HIDES at OAO and High Dispersion 
Spectrograph at Subaru telescope (Kambe et al. 2002). The cell is put in front of the 
entrance slit of the spectrograph and its temperature is controlled to 60 °C. We use a blue- 
arm and a middle-camera configuration which covers a wavelength range of 3900-7260A with 
a spectral resolution of ~ 40000 by two pixels sampling. Although a wide wavelength range 
can be obtained with a single exposure, only a waveband of 470 A is available for radial 
velocity measurements due to the small format (lKxlK) of the current CCD, which is a 
thinned back illuminated blue enhanced TEK CCD whose pixel size is 24x24 /xm 2 . 

The modeling technique of a star+I 2 spectrum and the extraction method of a stellar 
template are based on Sato et al. (2002) 's code. For CES data, we use 5 Gaussian profiles to 
reconstruct the instrumental profile and the second order Legendre polynomial function to 
express the wavelength scale. Entire echelle spectrum is divided into about forty segments 
(130-150 pixels for each), and Doppler analysis is applied to each segment. The current 
best short-term (~ 1 week) precision for bright stars is about 15 m s -1 limited by the low 
wavelength resolution and the narrow wavelength coverage of the spectrograph. In longer 
time scale (~ 1 yr), the precision gets worse to 20-25 m s -1 including velocity offsets between 
different observing runs. In the case of fainter stars, the measurement precision is mainly 
limited by a signal-to-noise ratio (S/N). We can typically achieve a precision of 30-40 m s _1 
for a V = 6 star with a S/N~150 for an exposure time of 1800 sec. 



Foundation, USA. 
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3. Stellar Properties 

11 Com (HR 4697, HD 107383, HIP 60202) is a G8 III giant star with a V magnitude 
V = 4.74, a color index B — V = 0.99, and a precise astrometric parallax it = 9.04 ± 0.86 
mas (ESA 1997). The resulting distance is 112 pc from the sun and the absolute magnitude 
is My = —0.48. Effective temperature T e g = 4742 ± 100K was derived from the B — V and 
metallicity [Fe/H] using the empirical calibration of Alonso et al. (2001), and a bolometric 
correction B.C. = —0.36 was derived from the calibration of Alonso et al. (1999) depending 
on temperature and metallicity. Then the stellar luminosity was estimated to L = 172±31L , 
and a radius to R = 19.5±2i? . Stellar mass M = 2.7±0.3M Q was estimated from the star's 
position in the theoretical H-R diagram by interpolating in the evolutionary tracks of Yonsei- 
Yale (Yi et al. 2003). Surface gravity log g = 2.31 ± 0.10 was determined using the relation 
between temperature, mass, luminosity and gravity (see equation 1 of Chen et al. 2000). Iron 
abundance was determined from the equivalent widths measured from the l2-free spectrum 
(5000-6100 and 6000-7100 A) combined with the model atmosphere (Kurucz 1993). We 
iterated the whole procedure described above until the value of the overall metallicity of 
the model converged. Finally we got [Fe/H] = —0.35 ± 0.09 and a microturbulent velocity 
Vt = 1.5 ± 0.2 km s -1 . The rotational velocity of the star is determined by de Medeiros & 
Mayor (1999) to be v sini = 1.2 ± 1.0 km s -1 . Hipparcos made a total of 99 observations of 
the star, revealing a photometric stability down to o = 0.006 mag. The stellar parameters 
are summarized in Table [TJ 

Figure 1 shows Ca II H line for 11 Com together with those for G8-9 III giants in our 
sample. Although the calibration to measure chromospheric activity for our spectra and the 
correlation between chromospheric activity and intrinsic radial velocity "jitter" for giants 
have not been well established yet, the lack of significant emission in the line of 11 Com 
suggests its low chromospheric activity. 

4. Radial Velocities and Orbital Solution 

We started observations of 11 Com at OAO in 2003 December. Soon after that, its large 
radial velocity variations were detected suggesting the existence of a substellar companion. 
Based on this result, we started observations of this star at Xinglong in 2005 to check the 
measurement precision with CES and to confirm the variations of the star independently. 
Until February 2007, we gathered a total of 28 data points of the star at OAO with a typical 
S/N of 250 per pixel, and a total of 18 data points at Xinglong with a typical S/N of 150-200. 

The observed radial velocities are shown in Figure [2] and listed in Table [2] (OAO) and 
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|3] (Xinglong) together with their estimated uncertainties. The best-fit Keplerian orbit was 
derived using both of the 0A0 and Xinglong data. An offset of —205 m s _1 was applied 
to the Xinglong radial velocities in order to minimize \ 2 when fitting a Keplerian model 
to the combined OAO and Xinglong velocities. The resulting orbit is shown in Figure [2] 
overplotted on the velocities, and its parameters are listed in Table HI The uncertainty of 
each parameter was estimated using a Monte Carlo approach. The radial velocity variability 
can be well fitted by a Keplerian orbit with a period P = 326.03 ± 0.32 days, a velocity 
semiamplitude K% = 302.8 ± 2.6 m s _1 , and an eccentricity e = 0.231 ± 0.005. When we 
determined a Keplerian orbit using only OAO data, the rms scatter of the residuals to the 
Keplerian fit for OAO data was 17.8 m s _1 , which is slightly larger than the typical intrinsic 
radial velocity scatter of late G-type giants in our sample (Sato et al. 2005), but we found no 
significant additional periodicity in the residuals at this stage. The Keplerian orbit derived 
by only OAO data fits the Xinglong data with the rms scatter of 33.2 m s _1 , which is 
comparable to the observational errors. Adopting a stellar mass of 2.7±0.3 M Q , we obtained 
for the companion a mass m2sini = 19.4 ± 1.5 Mj and a semimajor axis a = 1.29 ± 0.05 
AU. The uncertainties mostly come from that in host star's mass. If we assume the orbit 
is randomly oriented, there is a 3% chance that the true mass exceeds 80 Mj (i < 14°), a 
border between brown dwarf and star mass regimes. 

To investigate other causes producing apparent radial velocity variations such as pulsa- 
tion and rotational modulation, spectral line shape analysis was performed with use of high 
resolution stellar templates followed by the technique of Sato et al. (2007). We extracted 
two stellar templates from the star+I 2 spectra obtained at OAO: one was from five spectra 
with observed radial velocities of 250 ~ 300 m s _1 and the other was from those with ~ —200 
m s -1 . Cross correlation profiles of the templates were calculated for 50 spectral segments 
(4-5A width each) in which severely blended lines or broad lines were not included. Three 
bisector quantities were calculated for the cross correlation profile of each segment: the ve- 
locity span (BVS), which is the velocity difference between two flux levels of the bisector; 
the velocity curvature (BVC), which is the difference of the velocity span of the upper half 
and lower half of the bisector; and the velocity displacement (BVD), which is the average 
of the bisector at three different flux levels. We used flux levels of 25%, 50%, and 75% of 
the cross correlation profile to calculate the above quantities. Figure [3] shows the resulting 
bisector quantities plotted against the center wavelength of each segment. As expected from 
the planetary hypothesis, both of the bisector velocity span and the curvature are identical 
to zero (—2.3 m s _1 and —2.4 m s _1 in average, respectively), which means that the cross 
correlation profiles are symmetric, and the average bisector velocity displacement of —453.0 
m s _1 is consistent with the velocity difference between the two templates. Furthermore, 
these values are irrelevant to wavelength. Based on these results, we conclude that the radial 
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velocity variability observed in 11 Com is best explained by orbital motion. 

5. Discussion and Summary 

We reported a discovery of a brown dwarf mass companion to the intermediate-mass 
giant 11 Com. This is the first result from the joint planet search program between China 
and Japan. We obtained radial velocities of the star over a time span of 3 years at OAO 
and 2 years at Xinglong. The two sites independently detected the velocity variability of the 
star, and the agreement of the two data sets assures the reliability of our discovery and the 
radial velocity measurements at both sites. 

11 Com b is the third substellar companion ever discovered to intermediate- mass giants 
falling in the brown dwarf regime (FigHJ). To date, 10 planetary mass companions (0.6 — 
10Mj) have been found around intermediate-mass (> 1.6M ) giants by precise Doppler 
surveys, but only 3 brown dwarf mass ones (> 14Mj) have been found despite the larger 
radial velocity amplitudes being easily detectable. These results may suggest a shortage of 
brown dwarfs relative to planets around intermediate-mass giants within orbital separations 
of -3 AU. 

Brown dwarf companions are generally thought to form by gravitational collapse in 
protostellar clouds like stellar binary systems (Bonnell & Bastien 1992; Bate 2000). It is 
difficult, however, to form close binary systems with large difference in mass between the 
primary and the secondary (Bate 2000). According to this scenario, fewer brown dwarf 
companions are expected in higher-mass systems than in lower-mass ones. Around solar- 
type stars, the frequency of brown dwarf companions is estimated to be less than 0.5% 
within 3 AU (Marcy & Butler 2000) @. As long as our survey, 11 Com b emerged from a 
sample of 300 stars at OAO. The current detection rate seems to be comparable to that 
around solar-type stars, although the exact detection limits for all the targets in our sample 
have not been derived yet. Detailed statistical analysis will be presented in a forthcoming 
paper. Ongoing large scale Doppler surveys of G K giants (Setiawan et al. 2005; Sato et 
al. 2007; Hatzes et al. 2005; Johnson et al. 2007; Lovis & Mayor 2007; Frink et al. 2002; 
Niedzielski et al. 2007) and B-A dwarfs (Galland et al. 2005) around the globe will derive 
statistically reliable value of the frequency. 

Gravitational instability in protostellar disks is another scenario to form brown dwarf 



3 Patel et al. (2007) recently discovered five brown dwarf candidates beyond 3 AU by precise radial velocity 
technique. 
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companions (Boss 2000; Rice et al. 2003). This scenario is also proposed as the formation 
mechanism for massive gaseous planets. A half of the planetary mass companions discovered 
around intermediate-mass stars have minimum masses greater than 5 Mj, so-called "super 
planets" (Hatzes et al. 2005; Sato et al. 2003, 2007; Setiawan et al. 2003, 2005; Lovis & 
Mayor 2007). Can these companions and 11 Com b be put into the same context in the 
point of formation mechanism? One remarkable feature of these companions is that most 
of them (including 11 Com b) are orbiting stars with subsolar metallicity. This apparently 
favors gravitational instability scenario as their formation mechanism (Boss 2002), which is 
less sensitive to metallicity than core accretion scenario is (Fischer & Valenti 2005; Ida & 
Lin 2004). However, if the surface density of dust in the disk scales as Mf Z (a > 0), where 
M* is the stellar mass and Z is the metallicity, massive stars may have higher disk surface 
densities, which may compensate for lower metalicities. Such disks could form massive solid 
cores that can accrete huge gas envelope up to 20Mj. Further investigation of mass and 
metallicity distribution of substellar companions to massive stars helps discriminate between 
these formation scenarios. 

This research is based on data collected at Xinglong station, which is operated by the 
National Astronomical Observatories, CAS (NAOC), and Okayama Astrophysical Observa- 
tory (OAO), which is operated by National Astronomical Observatory of Japan (NAOJ). 
We thank Junjun Jia, Jiaming Ai, and Hongbin Li for their expertise and support of the 
Xinglong observations. We are grateful to all of the staffs of OAO for their support during 
the OAO observations. Data analysis was in part carried out on "sb" computer system op- 
erated by the Astronomical Data Analysis Center (ADAC) and Subaru Telescope of NAOJ. 
We thank the National Institute of Information and Communications Technology for their 
support on high-speed network connection for data transfer and analysis. This work was 
funded by the National Natural Science Foundation of China under grants 10521001 and 
10433010. B.S., H.I., H.A., and M.Y. are supported by Grants-in-Aid for Scientific Research 
No.17740106, (C) No.13640247, (B) No.17340056, (B) No.18340055, respectively, from the 
Japan Society for the Promotion of Science (JSPS). E.T., D.M., and Y.I. are supported 
by "The 21st Century COE Program: The Origin and Evolution of Planetary Systems" in 
Ministry of Education, Culture, Sports, Science and Technology (MEXT). E.K. and S.I. are 
partially supported by MEXT, Japan, the Grant-in-Aid for Scientific Research on Priority 
Areas, "Development of Extra-Solar Planeary Science." This research has made use of the 
SIMBAD database, operated at CDS, Strasbourg, France. 
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Fig. 1. — Spectra in the region of Call H lines. HD 120048 and HD 176598 show significant 
core reversals in the lines suggesting high chromospheric activity. They have relatively large 
radial velocity scatters of 25-35 m s^ 1 . 7 Tau exhibits no emission in the line and show 
radial velocity stability at a level of ~5 m s _1 . 
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Fig. 2. — Radial velocities of 11 Com (dots) observed at OAO (red) and Xinglong (green). 
The Keplerian orbit (solid line) is determined using both of the OAO and Xinglong data. 
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Fig. 3. — Bisector quantities of cross-correlation profiles between the templates of 11 Com at 
peak (250-300 m s _1 ) and valley (—200 m s" 1 ) phases of observed radial velocities, showing 
bisector velocity span (BVS, circle), bisector velocity curvature (BVC, triangle), and bisector 
velocity displacement (BVD, square). The definition of these quantities are described in 
Section 4. The mean values and their standard errors are shown in the figure. Dashed lines 
represent mean values of the quantities (an offset of 300 m s _1 is added to the BVS). 
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Fig. 4. — Mass distribution of substellar companions discovered by precise Doppler sur- 
veys (including all of the host stars with < 4M Q ). The data are based on the list from 
|http://exoplanets.org/planets.shtml| NGC2423 No3 b (10.6 Mj), NGC4349 Nol27 b (19.8 
Mj), HD 13189 b (K2 II, 14Mj), HD 17092 b (K0 III, 4.6 Mj) and HD 180777 b (A9 V, 
25Mj) are not included in the list but are added to the histogram. 
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Table 1. Stellar Parameters for 11 Com 



Parameter 


Value 


Sp. Type 


G8 III 


7r (mas) 


9.04±0.86 


V 


4.74 


B-V 


0.99 


My 


-0.48 


B.C. 


-0.36 


T cS (K) 


4742±100 




2.31±0.1 


[Fe/H] 


-0.35±0.09 


Vt (km s _1 ) 


1.5±0.2 


L(L & ) 


175±31 


R (R e ) 


19±2 


M (M ) 


2.7±0.3 


v sini (km s _1 ) 


1.2±1.0 
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Table 2. Observed Radial Velocities of 11 Com at OAO 



JD 


Radial Velocity 


Error 


(-2,450,000) 


(m s _1 ) 


(m s _1 ) 


3002.2791 


-201.7 


5.7 


3101.0241 


178.7 


4.6 


3216.0035 


74.5 


7.9 


3334.3172 


-184.5 


5.5 


3340.3401 


-160.2 


4.6 


3364.3191 


-76.9 


7.9 


3366.2249 


-76.6 


5.2 


3367.2333 


-79.3 


5.4 


3376.3243 


-12.5 


4.7 


3403.2636 


65.3 


5.1 


3424.1132 


150.7 


5.1 


3447.1803 


232.5 


4.6 


3468.0504 


251.2 


5.6 


3495.1120 


293.4 


4.5 


3521.0884 


212.5 


7.5 


3577.9651 


-217.7 


6.2 


3694.3575 


-38.2 


5.0 


3728.3204 


63.2 


5.1 


3743.3484 


127.7 


5.0 


3775.2089 


241.6 


4.7 


3810.1629 


311.6 


3.9 


3831.1262 


301.0 


4.5 


3853.1393 


188.7 


5.8 


3887.0520 


-148.5 


5.7 


4051.3576 


43.8 


4.6 


4088.3683 


189.0 


4.8 


4121.2706 


239.6 


5.9 


4143.2374 


274.5 


4.0 
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Table 3. Observed Radial Velocities of 11 Com at Xinglong 



JD 


Radial Velocity 


Error 


(-2,450,000) 


(m s _1 ) 


(m s _1 ) 


3420.2939 


195.4 


33.5 


3422.2517 


126.7 


36.2 


3422.2627 


112.0 


33.1 


3511.0423 


281.9 


29.9 


3514.0647 


249.4 


25.6 


3514.0779 


304.8 


20.6 


3775.2678 


243.1 


30.6 


3778.3433 


220.6 


30.3 


3839.1013 


229.1 


36.6 


3839.1229 


204.4 


18.6 


3839.1462 


246.7 


31.8 


3841.1569 


227.5 


28.9 


3891.0521 


-124.4 


18.3 


3894.0998 


-178.1 


22.7 


4045.3959 


78.7 


21.1 


4046.3706 


58.3 


26.7 


4132.3744 


249.1 


30.9 


4132.3976 


253.2 


19.7 
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Table 4. Orbital Parameters of 11 Com determined from both of 0A0 and Xinglong 



data. 


Parameter 


Value 


P (days) 


326.03±0.32 


K x (m s- 1 ) 


302.8±2.6 


e 


0.231±0.005 


u (deg) 


94.8±1.5 


T p (JD-2,450,000) 


2899.6±1.6 


aisini (10~ 3 AU) 


8.848±0.073 


/i(m) (1(T 7 M ) 


8.68±0.21 


Tasini (Mj) 


19.44:1.5 


a (AU) 


1.29±0.05 


N 


46 


rms (m s _1 ) 


25.5 


Reduced \fy? 


2.8 



